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Abstract 

In this paper, we consider secure downlink transmission in a multi-cell massive multiple-input 
multiple-output (MIMO) system where the numbers of base station (BS) antennas, mobile terminals, 
and eavesdropper antennas are asymptotically large. The channel state information of the eavesdrop¬ 
per is assumed to be unavailable at the BS and hence, linear precoding of data and artificial noise 
(AN) are employed for secrecy enhancement. Four different data precoders (i.e., selfish zero-forcing 
(ZF)/regularized channel inversion (RCI) and collaborative ZF/RCI precoders) and three different AN 
precoders (i.e., random, selfish/collaborative null-space based precoders) are investigated and the cor¬ 
responding achievable ergodic secrecy rates are analyzed. Our analysis includes the effects of uplink 
channel estimation, pilot contamination, multi-cell interference, and path-loss. Furthermore, to strike a 
balance between complexity and performance, linear precoders that are based on matrix polynomials are 
proposed for both data and AN precoding. The polynomial coefficients of the data and AN precoders 
are optimized respectively for minimization of the sum mean squared error of and the AN leakage to the 
mobile terminals in the cell of interest using tools from free probability and random matrix theory. Our 
analytical and simulation results provide interesting insights for the design of secure multi-cell massive 
MIMO systems and reveal that the proposed polynomial data and AN precoders closely approach the 
performance of selfish RCI data and null-space based AN precoders, respectively. 


I. Introduction 

Massive multiple-input multiple-output (MIMO) systems employing simple linear precoding 
and combining schemes offer significant performance gains in terms of bandwidth, power, and 
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energy efficiency compared to conventional multiuser MIMO systems as impairments such as 
fading, noise, and interference are averaged out for very large numbers of base station (BS) 
antennas m-m. Furthermore, in time-division duplex (TDD) systems, channel reciprocity can 
be exploited to estimate the downlink channels via uplink training so that the training overhead 
scales only linearly with the number of users and is independent of the number of BS antennas 
[[2l • However, if the pilot sequences employed in different cells are not orthogonal, so-called pilot 
contamination occurs and impairs the channel estimates, which ultimately limits the achievable 
performance of massive MIMO systems l(2ll. ||4||. 

Since secrecy and privacy are critical concerns for the design of future communication systems 
|[5l , it is of interest to investigate how the large number of spatial degrees of freedom in 
massive MIMO systems can be exploited for secrecy enhancement |6]|, 0. If the eavesdropper 
(Eve) remains passive to hide its existence, neither the transmitter (Alice) nor the legitimate 
receiver (Bob) will be able to leam Eve’s channel state information (CSI). In this situation, it 
is advantageous to inject artificial noise (AN) at the transmitter to degrade Eve’s channel and 
to use linear precoding to avoid impairment to Bob’s channel as was shown in (8l- lITOll and 
ifTTl . llT2l for single user and single-cell multiuser systems, respectively. However, in multi-cell 
massive MIMO systems, multi-cell interference and pilot contamination will hamper Alice’s 
ability to degrade Eve’s channel and to protect Bob’s channel. This problem was studied first 
in lfl3l for simple matched-filter (MF) data precoding and null-space (NS) and random AN 
precoding. However, it is well known that MF data precoding suffers from a large loss in the 
achievable information rate compared to other linear data precoders such as zero-forcing (ZF) and 
regularized channel inversion (RCI) precoders as the number of mobile terminals (MTs) increases 
DU- Since it is expected that this loss in information rate also translates into a loss in secrecy 
rate, studying the secrecy performance of ZF and RCI data precoders in massive MIMO systems 
is of interest. Furthermore, while NS AN precoding was shown to achieve a better performance 
compared to random AN precoding llT3ll . it also entails a much higher complexity. Similarly, the 
improved performance of ZF and RCI data precoding compared to MF data precoding comes at 
the expense of a higher complexity. Hence, the design of novel data and AN precoders which 
allow a flexible tradeoff between complexity and secrecy performance is desirable. 

Related work on physical layer security in massive MIMO systems includes llT5l where the 
authors use the channel between Alice and Bob as secrete key and show that the complexity 
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required by Eve to decode Alice’s message is at least of the same order as a worst-case lattice 
problem. Physical layer security in a downlink multi-cell MIMO system was considered in lfT6l - 
|fT8Tl . However, unlike our work, perfect knowledge of Eve’s channel was assumed, AN injection 
was not considered, and pilot contamination was not taken into account. Furthermore, ZF and 
RCI data precoding were analyzed in the large system limit in |fl9l , l!20l . However, neither pilot 
contamination nor AN were taken into account and the secrecy rate was not analyzed. Using a 
concept that was originally conceived for code division multiple access (CDMA) uplink systems 
in ll2Tll and later extended to MIMO systems in ll22l reduced complexity linear data precoders 
that are based on matrix polynomials were investigated for use in massive MIMO systems in 
lf23l - ll25ll . However, ll23l - lf25ll did not take into account the effect of AN leakage for precoder 
design and did not study the secrecy performance. Hence, the results presented in IT5l - ll25ll are 
not directly applicable to the system studied in this paper. 

In this paper, we consider secure downlink transmission in a multi-cell massive MIMO system 
employing linear data and AN precoding in the presence of a passive multi-antenna eavesdropper. 
We study the achievable ergodic secrecy rate of such systems for different linear precoding 
schemes taking into account the effects of uplink channel estimation, pilot contamination, multi¬ 
cell interference, and path-loss. The main contributions of this paper are summarized as follows: 

• We study the performance-complexity tradeoff of selfish and collaborative data and AN 
precoders. Selfish precoders require only the CSI of the MTs in the local cell but cause 
inter-cell interference and inter-cell AN leakage. In contrast, collaborative precoders require 
the CSI between the local BS and the MTs in all cells, but reduce inter-cell interference 
and inter-cell AN leakage. However, since the additional CSI required for the collaborative 
precoders can be estimated directly by the local BS, the additional overhead and complexity 
incurred compared to selfish precoders is limited. 

• We derive novel closed-form expressions for the asymptotic ergodic secrecy rate which 
facilitate the performance comparison of different combinations of linear data precoders 
(i.e., MF, selfish and collaborative ZF/RCI) and AN precoders (i.e., random, selfish and 
collaborative NS), and provide significant insight for system design and optimization. 

• In order to avoid the computational complexity and potential stability issues in fixed point 
implementations entailed by the large-scale matrix inversions required for ZF and RCI data 
precoding and NS AN precoding, we propose polynomial (POFY) data and AN precoders 
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and optimize their coefficients. Unlike Il24ll and 11251 . which considered polynomial data 
precoders for massive MIMO systems without AN generation, we use free probability theory 
Il23l , ||26l to obtain the POLY coefficients. This allows us to express the POLY coefficients 
as simple functions of the channel and system parameters. Simulation results reveal that 
these precoders are able to closely approach the performance of selfish RCI data and NS 
AN precoders, respectively. 

The remainder of this paper is organized as follows. In Section [TT1 we outline the considered 
system model and review some basic results from lTT3l . In Sections HU and |W] the considered 
linear data and AN precoders are investigated, respectively. In Section [V] the ergodic secrecy rates 
of different linear precoders are compared analytically for a simple path-loss model. Simulation 
and numerical results are presented in Section [VH and some conclusions are drawn in Section 

m 

Notation: Superscripts T and H stand for the transpose and conjugate transpose, respectively. 
Ijv is the A-dimensional identity matrix. The expectation operation and the variance of a random 
variable are denoted by E[-] and var[-], respectively. diag{x} denotes a diagonal matrix with the 
elements of vector x on the main diagonal. tr{ -} and rank{-} denote trace and rank of a matrix, 
respectively. C mxri represents the space of all m x n matrices with complex-valued elements, 
x ~ CN(Oat, S) denotes a circularly symmetric complex Gaussian vector x e C /Vx 1 with zero 
mean and covariance matrix £. [A]^ denotes the element in the k th row and I th column of 
matrix A, and [x] + = max{x, 0}. 

II. System Model and Preliminaries 

In this section, we introduce the considered system model as well as the adopted channel 
estimation scheme, and review some ergodic secrecy rate results. 

A. System Model 

We consider the downlink of a multi-cell massive MIMO system with M cells and a frequency 
reuse factor of one, i.e., all BSs use the same spectrum. Each cell includes one Ay-antenna BS, 
K < Ay single-antenna MTs, and potentially an Ay-antenna eavesdropper. The eavesdroppers 
try to hide their existence and hence remain passive. As a result, the BSs cannot estimate the 
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eavesdroppers’ CSI. To overcome this limitation, each BS generates AN to mask its information¬ 
carrying signal and to prevent eavesdropping [[§]]. In the following, the k th MT, k — 1,... ,K, 
in the n th cell, n = 1,..., M, is the MT of interest and we assume that an eavesdropper tries to 
decode the signal intended for this MT. We note that neither the BSs nor the MTs are assumed 
to know which MT is targeted by the eavesdropper. The signal vector, x n G C /VrX 1 , transmitted 
by the BS in the n th cell (also referred to as the n ih BS in the following) is given by 

-X n \/pF n S n T ^/f^A n Z n , (1) 

where s„ ~ CN(0 ,a, I a) and z„ ~ CN(0 at t . I \ 7 .) denote the data and AN vectors for the K MTs 
in the n th cell, respectively. F„ = [f nl , ■ ■ • , f nK ] G C NtxK and A n = [a nl , • • • , a nNr } G C NtxNt 

are the data and AN precoding matrices, respectively, and the efficient design of these matrices 
is the main scope of this paper. Thereby, the structure of both types of precoding matrices does 
not depend on which MT is targeted by the eavesdropper. The AN precoding matrix A n has 
rank L = rank{ A n } < N E , i.e., L dimensions of the Ag-dimensional signal space spanned by 
the Nr BS antennas are exploited for jamming of the eavesdropper. The data and AN precoding 
matrices are normalized as tr{F^F,„} = K and tr{A^A n } = L, i.e., their average power per 
dimension is one. The average powers p and cj allocated to the information-carrying signal for 
each MT and each AN signal, respectively, can be written as p = and q = -—where 
Pt is the total transmit power and o G (0,1] is a power allocation factor which can be optimized. 
For the sake of clarity, in this paper, we assume that all cells utilize the same value of cj). 

The vectors collecting the received signals at the K MTs and the Ne antennas of the 
eavesdropper in the n 1h cell are given by 

M M 

y n = ^T G mn x m + n n and y E = ^ G rnE x m + n E , (2) 

m =1 m— 1 

respectively, with Gaussian noise vectors n„ e CKf(0/<, and n E G CN(0^ £ , a|I jVfi ), 

where <7„ an ^ a E denote the noise variances at one MT and one eavesdropper receive antenna, 
respectively. Furthermore, G mn = D™H mn G C KxNt and G mE = G C NexNt are 

the matrices modeling the channels from the m th BS to the K MTs and the eavesdropper in 
the n th cell, respectively. Thereby, D mri = cliag{/3^ n ,.... 6^ n } and l3 rnE represent the path- 
losses from the m th BS to the K MTs and the eavesdropper in the n ih cell, respectively. Matrix 
H mn G C KxNt , with row vector h^ n G C 1xNt in the k th row, and matrix H mE G C NexNt 
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represent the corresponding small-scale fading components. Their elements are modeled as 
mutually independent and identically distributed (i.i.d.) complex Gaussian random variables 
(RYs) with zero mean and unit variance. 

For the design of the data and noise precoders, we consider two different approaches: Selfish 
designs and collaborative designs. For the selfish designs, each BS designs its precoders only 
based on the estimate of the CSI in its own cell, G nn , and without regard for the interference 
and the AN it causes to other cells. In contrast, for the collaborative designs, each BS designs 
its precoders based on the estimates of the CSI to the MTs in all cells, G mn , m — 1,..., M, in 
an effort to avoid excessive interference and AN to other cells. Although collaborative designs 
introduce more channel estimation overhead at the BS, they may not always outperform selfish 
designs because of the imperfection of the CSI and the limited number of spatial degrees of 
freedom available for precoder design. 


B. Channel Estimation and Pilot Contamination 


As is customary for massive MIMO systems, we assume that the downlink and uplink channels 
are reciprocal and the CSI is estimated in an uplink training phase Q]|- SI- To this end, all MTs 
emit pilot sequences of length r > K and with pilot symbol power p T . We assume that the pilot 
sequences of the K MTs in a given cell are mutually orthogonal but the same pilot sequences 
are used in all cells. This gives rise to so-called pilot contamination [0- 0. Furthermore, we 
assume that the path-loss information changes on a much slower time scale than the small-scale 
fading. Hence, the path-loss matrices D nm , m — 1,..., M, can be estimated perfectly and are 
assumed to be known at the BS for minimum mean-square error (MMSE) estimation of the 
small-scale fading gains 0. At the n' h BS, the small-scale fading vector to the k ih MT in the 
m th cell, hf' m , can be expressed as 


h 


k 

nm 


h nm. + h 


k 

nm'i 


(3) 


where the estimate lij'; m and the estimation error h k nm are mutually independent and can be statisti¬ 


cally characterized as h|; 


CN(0 


Nt 


PtTP 




Nt. 


and h 1 ' 


1 +P- rT S2l^m PnlT \ 

LN(0jVr ’ 1+P.TE&& ^ 


respectively, cf. Ifl3l . For future reference, we collect the estimates and the estimation errors at the 
BS corresponding to all K MTs in the m th cell in matrices H nm = [(h^ m ) T ,..., (h^ m ) T ] r £ 
C KxNt and H nm = [(X m ) T ,..., (h% m ) T ] T e C KxNt , respectively. 


th 
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C. Ergodic Secrecy Rate 

The performance metric adopted in this paper is the ergodic secrecy rate Q. In this section, we 
review some results for the ergodic secrecy rate in multi-cell massive MIMO systems employing 
linear data and AN precoding from lfT3ll . as these results will be needed throughout this paper. 
Combining © and ([2]) we observe that the downlink channel comprising the BS, the k th MT, 
and the eavesdropper in the n th 

cell is an instance of a multiple-input, single-output, multi¬ 
eavesdropper (MISOME) wiretap channel [[6]|. Hence, the achievable secrecy rate of the k th MT 
in the n th cell is bounded by the difference of the capacities of the channel between the BS and 
the MT and the channel between the BS and the eavesdropper, see lfl3l Lemma 1], lfl7l Lemma 
2]. Thus, a lower bound on the ergodic secrecy rate of the k th MT in the n th cell is given by 

m 

R% = {R, k -CZ] + ,k=l,...,K, (4) 


where R n k denotes an achievable rate of the k th MT in the n th cell and C®™ denotes the ergodic 
capacity of the channel between the BS and the eavesdropper. In order to obtain a tractable 
lower bound on the ergodic secrecy rate, we lower bound the achievable rate of the MT as 
Rnk = log 2 (l + Ink) with signal-to-interference-and-noise ratio (SINR) [fT3l Eq. (10)] 

„ _ _ MVPnnP h nJnk}\ 2 _ 

' nk _ M N t 

™4VPLP h nJnk] + EE E (l\//S h mn a ™| 2 ] + E E [| V ^nP h tJml\ 2 } + 1 

m= 1 i —1 {m,l}^{n,k} 

(5) 


Lurthermore, we make the pessimistic assumption that the eavesdropper is able to cancel the 
received signals of all in-cell and out-of-cell MTs except the signal intended for the MT of 
interest. This leads to an upper bound for the eavesdropper’s capacity, and consequently, to a 
lower bound for the ergodic secrecy rateQ Hence, the ergodic capacity of the eavesdropper is 
given by dll Eq. (7)] 


C eve 
nk 


= E 


log; 


(l+p&G^X- 


' G n E^nk) 


( 6 ) 


where X = ?Em=i G m £A m A^G^ E e £N t xn t ^ enotes the noise correlation matrix at the 
eavesdropper under the worst-case assumption that the receiver noise at the eavesdropper is 


'This lower bound is achievable if the eavesdropper has access to the data of all interfering in-cell and out-of-cell MTs, which 
might be the case e.g. if the interfering MTs cooperate with the eavesdropper. 
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negligible, i.e., a% —> 0. Denoting the normalized number of eavesdropper antennas by a = 
N e /N t , a necessary condition for the invertibility of matrix X is a < ML/N T . Hence, a non¬ 
zero secrecy rate can only be achieved if this condition is met. Consequently, a larger L implies 
that the BS is able to tolerate more eavesdropper antennas. 

If H nS f„ fc and matrix X are statistically independent, which in turn means for the data and 
AN precoders that vector f n k and the subspace spanned by the columns of A„ are mutually 
orthogonal, a simple and tight upper bound on © can be obtained. Since any efficient data/AN 
precoder pair has to keep the AN self-interference at the desired MT small, this orthogonality 
condition holds at least approximately in practice. In this case, for a < a 2 L/(cN T ) and N T —> 
oo, where a = 1 + P™e / PnE and c = 1 + / Pue) 2 , a simple and tight upper 

bound for C° v k e is given by lfl3l Theorem 1] 

° nk < !og 2 (f + aqL / Nrr _ caq /^j = log 2 ( X + p(i - 0)(a - caN T /(La)y) ' (7) 

For M — 1, we have a 2 /c = M = 1, i.e., the bound in (0) is applicable in the entire 
range of a where C^ in © is finite. For M > 1, we have a 2 /c < M, i.e., the bound is not 
applicable for La 2 /(cN T ) < a < ML/N T . However, for strong inter-cell interference, we have 
/3 m E ~ /3 u e and a 2 /c ~ M, i.e., the bound is applicable for all a for which C^ k e in © is 
finite. On the other hand, for weak inter-cell interference, we have /3 m E <C (3 n E, and matrix 
X will be ill-conditioned for L/N T < a < ML/N T and will become very large. Hence, 
the bound is again applicable for the values of a (i.e., 0 < a < L/N T ), for which C®™ in 
© assumes practically relevant values. More generally, [T3\ Figs. 2-4] and Section [Vi] suggest 
that, for N t oc. (0) is applicable and tight for all values of a which permit a non-vanishing 
secrecy rate. 

Combining ©, ©, and 0). we obtain a tight and tractable lower bound on the secrecy rate 
llT3l . It is noteworthy that the upper bound on the capacity of the eavesdropper in © is only 
affected by the dimensionality of the AN precoder, L, but not by the exact structures of A n and 
F n , as long as i n k and the subspace spanned by the columns of A n are orthogonal. On the other 
hand, the achievable rate of the MT in © is affected by both the data and the AN precoders. 
In the following two sections, we analyze the impact of the most important existing data and 
AN precoder designs on the achievable rate R n k as N T -A oo, respectively, and propose novel 
low-complexity data and AN precoders that are based on a polynomial matrix expansion. 
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III. Linear Data Precoders for Secure Massive MIMO 


In this section, we analyze the achievable rate of selfish and collaborative ZF/RCI data 
precoding, respectively, and develop a novel POLY data precoder. In contrast to existing analyses 
and designs of data precoders for massive MIMO, e.g. Ifl9l . Il20ll . Il23ll - ll25ll . the results presented 
in this section account for the effect of AN leakage, which is only present if AN is injected at 
the BS for secrecy enhancement. We are interested in the asymptotic regime where K , Nr —» oo 
but 3 = K/N t and a = N E /N T are finite. 


A. Analysis of Existing Data Precoders 

For Nt —> oo, analyzing the achievable rate is equivalent to analyzing the SINR in ©. 
Thereby, the effect of the AN precoder can be captured by the term 

M Nt M 

Q = E E E[l v&La™! 2 ] = E A„E[hLA,,.A"(hL) H ] (8) 

m= 1 2=1 771=1 

in the denominator of ©, which represents the inter-cell and intra-cell AN leakage. This term 
is assumed to be given in this section and will be analyzed in detail for different AN precoders 
in Section [TV] 

1) Selfish ZF/RCI Data Precoding: The selfish RCI (SRCI) data precoder for the n lh cell is 
given by 

F„ = 7, L nn H" n , (9) 


where L nn = (H^ n H ran + n i I y r ) ~ 1 , 71 is a scalar normalization constant, and K\ is a regular¬ 
ization constant. In the following proposition, we provide the resulting SINR of the k th MT in 
the n th cell. 

Proposition 1 : For SRCI data precoding, the received SINR at the k th MT in the n th cell is 
given by 


7nfc CI 


f S RCI + (l+g(/9Al)) 2 


g(^i) fsRCi+ rgR f K1 (l+g(/3Ai)) 2 


4- W B k IB 

1 / jrriy^n ^mn/ 


3 k 

nn 


( 10 ) 


where 


= 2 


( 1 ~ ft) 2 , 2(1+ J9) | 1 | 1-/3 




Ki 


K1 


(ID 
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SRCI 


E m ^ 

m^n Z_/ 


l^k / 


fj , — _ PT T {Pmn) _ ,,Q — Pk 


+^+^’ i+p^Er=iC, 


and f SR0I = r gXli wi(h r ~— “-— 

s, 1 + P' rT E^ m /3L , _ / 

X l+PrrEH 1 ^ l ’ and ^-^- 

Proof: Please refer to Appendix A. ■ 

Regularization constant K\ can be optimized for maximization of the lower bound on the 

secrecy rate in ©, which is equivalent to maximizing the SINR in (ITOb . Setting the derivative 

of 7 ®fc CI with respect to k\ to zero, the optimal regularization parameter is found as k \, opt = 

/3/r SRC i, and the corresponding maximum SINR is given by 

1 


7nfc CI = 


( 12 ) 


l/^(/3, K l,opt) + fimn/finn 

On the other hand, for K\ —> 0, the SRCI data precoder in © reduces to the selfish ZF (SZF) 
data precoder. The corresponding received SINR is provided in the following corollary. 

Corollary 1 : Assuming /3 < 1, for SZF data precoding, the received SINR at the k ih MT in 
the n th cell is given by 

(13) 


7 nfc F = 


- 4 - V B k / B k 

1 / jm^n r^mn/ r^nn 


V rySRCI 


(1—/3)fsRCi ' ^rn^n > 

Proof: in (fl3l) can be obtained from (fTOl) as 7 ®f F = lim* 

2) Collaborative ZF/RCI Precoding: The collaborative RCI (CRCI) precoder for the n th cell 
is given by 

F„ = 72L„HL (14) 


where L n = (H^H n + k 2 Iat t ) _1 with H„ = [H F X ... H F M ] 7 G C MKxNt , 72 is a normalization 
constant, and k 2 is a regularization constant. The corresponding SINR of the k th MT in the n lh 
cell is provided in the following proposition. 

Proposition 2: For CRCI data precoding, the received SINR at the k th MT in the n th cell is 
given by 


7nfc RCI 


r C Rci+(i+g(M/3, K2 ))2 


g(M/3,* 2 ) ( r C RCi+ rcR ? lK2 (1+g(M/3, K2 )) 2 




k 

nn 


(15) 


where f CR ci = fSSfel with r ™ci = JfA 

Proof: The proof is similar to that for the SINR for the SRCI data precoder given in 
Appendix A and omitted here for brevity. ■ 
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Furthermore, the optimal regularization constant maximizing the SINR (and thus the secrecy 
rate) in (fl5T) is obtained as k 2 ,opt = M/3 /Tcrci, and the corresponding maximum SINR is given 
by 


7„T CI = 


( 16 ) 


1/G{MP, K 2 ,opt) + 

On the other hand, for k 2 —> 0, the CRCI precoder in (fl4l) reduces to the collaborative ZF 
(CZF) precoder. The corresponding received SINR is provided in the following corollary. 

Corollary 2: Assuming /? < 1/M, for CZF data precoding, the received SINR at the k th MT 
in the n th cell is given by 


Ink F = 


M/3 


(l—M^)r C Rci 


+ Pmn /A 


Jfc 

nn 


(17) 


Proof: 7nfc F in (IT7l) is obtained by letting k 2 — > 0 in dT31) . ■ 

Remark 1: Selfish data precoders require estimation of in-cell CSI, i.e., H nn , only. In contrast, 
collaborative data precoders require estimation of both in-cell and inter-cell CSI at the BS, i.e., 
H n . Furthermore, since collaborative data precoders attempt to avoid interference not only to in¬ 
cell users but also to out-of-cell users, more BS antennas are needed to achieve high performance. 
This is evident from Corollaries 1 and 2, which reveal that N T > K and N T > MK are 
necessary for SZF and CZF data precoding, respectively. On the other hand, if successful, trying 
to avoid out-of-cell interference is beneficial for the overall performance. Hence, whether selfish 
or collaborative precoders are preferable depends on the parameters of the considered system, 
cf. Sections [V] and [VU 


B. Polynomial Data Precoder 

The RCI and ZF data precoders introduced in the previous section achieve a higher perfor¬ 
mance than simple MF data precoding lfl3ll . However, they require a matrix inversion which 
entails a high computational complexity for the large values of K and Nr desired in massive 
MIMO. Hence, in this section, we propose a low-complexity POLY data precoder which avoids 
the matrix inversion. As the goal is a low-complexity design, we focus on selfish POLY precoders, 
although the extension to collaborative designs is possible. 
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The proposed POLY precoder, F n , for the n th BS can be expressed as 

x 


F = 

n x/iVr 


1 V// 

n tm / , Ti 


«=0 


_~_H 

H nn H nn 


(18) 


where H nn = ~^=H nn , and /i = [/i 0 ,.... //i ] 7 are the real-valued coefficients of the precoder 
matrix polynomial, which have to be optimized. In the following, we show that, for K, Nt —>■ oo, 
the optimum coefficients /i do not depend on the instantaneous channel estimates but are constant 
and can be determined by exploiting results from free probability lf26ll and random matrix theory 
ll29l . To this end, we define the asymptotic average mean-square error (MSE) of the users in the 
n th cell as mse n = lim k^oo 77 E [||e„|| 2 ] with error vector 


n c(G nn (-y/pF n S n T ^/^A n z n ) T n n | s n , (19) 

where n„ = Yhm^n G mn x m + n„ includes Gaussian noise, inter-cell interference, and inter-cell 
AN leakage. Furthermore, c is a normalization constant at the receiver, which does not impact 
detection performance. The optimal coefficient vector /x minimizes mse n for a given power 
budget 4>Pt for the information-carrying signal, i.e., 


min Mi? mse n s.t. : Tr{Ff F n } = 1, 


( 20 ) 


where we use the notation Tr {•} = lim^-Kx) tr {•} / K. The optimal coefficient vector, /x t , is 
provided in the following theorem. 

Theorem 1\ For K, N T —> oo, the optimal coefficient vector minimizing the asymptotic average 
MSE of the users in the n 1h cell for the POFY precoder in (fl8l) is given by 


/Apt = 73II 1 ' 0 , 


( 21 ) 


where -0 = [C, C 2 , - - -, C X+1 ] T , My = Tr {D„„} C i+j + (Tr {D nn A n } + ^ {S ^+ Pan ) C +j ~\ 
= E[n n nf], A n = diag 


1 +p T r J2m=l Pnn 


’ 1+PrT 


K 

nm 


NtP J 

and Pan = 


qE [Tr {G nn A n A^G^}]. Furthermore, ( l denotes the / th -order moment of the sum of the 

^ /s ff 

eigenvalues of H nn H nn , i.e., ( l = lim^oo £ Ef=i which converges to (' = Y!i =J Q (0J f 
for A" —)■ 00 ll23l Theorem 2]. Finally, 73 is chosen such that Tr{F^F n } = 1 holds. 

Proof: Please refer to Appendix B. ■ 

We note that /x opt does not depend on instantaneous channel estimates, and hence, can be 
computed offline. 
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C. Computational Complexity of Data Precoding 

We compare the computational complexity of the considered data precoders in terms of the 
number of floating point operations (FLOPs) lf28l . Each FLOP represents one scalar complex 
addition or multiplication. We assume that the coherence time of the channel is T symbol intervals 
of which r are used for training and T — r are used for data transmission. Hence, the complexity 
required for precoding in one coherence interval is comprised of the complexity required for 
generating one precoding matrix and T — r precoded vectors. A similar complexity analysis was 
conducted in (23l Section IV] for selfish data precoders without AN injection at the BS. Since 
the AN injection does not affect the structure of the data precoders, we can directly adapt the 
results from (23l Section IV] to the case at hand. In particular, the selfish MF, the SZF/SRCI, 
and the CZF/CRCI precoders require (2 K — 1 )N t (T — r), 0.5(A' 2 + K)(2N T — 1) + K 3 + K 2 + 
K + N t K(2K -1) + (2AT - 1 )N t (T - r), and 0.5(M 2 A' 2 + MAT) (2 N T -1) + M 3 K 3 + M 2 K 2 + 
MK + NtMK(2MK — 1) + (2AT — Y)Nj(T — r) FLOPs per coherence interval, see (231 Section 
IV]. In contrast, for the POLY data precoder, we obtain for the overall computational complexity 
(T — r) ((X + 1)(2A" — 1 )N t + X(2 N t — 1)AT) FLOPs, which assumes implementation of the 
precoding operation by Horner’s rule (231 Section IV]. 

The above complexity expressions reveal that the additional complexity introduced by collab¬ 
orative data precoders compared to selfish data precoders is at most a factor of M 3 . In addition, 
the complexity savings achieved with the POLY data precoder compared to the SZF/SRCI data 
precoders increase with increasing K for a given T. We note however that, regardless of their 
complexity, POLY data precoders are attractive as they avoid the stability issues that may arise 
in fixed point implementation of large matrix inverses. 

IV. Linear AN Precoders for Secure Massive MIMO 

In this section, we investigate the performance of selfish and collaborative NS (S/CNS) and 
random AN precoders. In addition, a novel POLY AN precoder is derived. To the best of the 
authors’ knowledge, POLY AN precoding has not been considered in the literature before. 

A. Analysis of Existing AN Precoders 

For a given dimensionality of the AN precoder, L, the secrecy rate depends on the AN precoder 
only via the AN leakage, Q, given in ([X]). which affects the SINR of the MT. Furthermore, the 
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optimal POLY data precoder coefficients in (I2T1) are affected by the AN precoder via the leakage 
term Pan- In this subsection, for N T —>■ oo, we will provide closed-form expressions for 0 and 
Pan for the SNS, CNS, and random AN precoders. 

1) SNS AN Precoder: The SNS AN precoder of the n th BS is given by [j8] 

-l „ 


A,„ = I 


n t 


H" ( U nn U nn 


f H 


Hr 


( 22 ) 


which has rank L = Nt — K and exists only if 13 < 1. We divide the corresponding AN 
leakage Qsns into an inter-cell AN leakage Q^ NS and an intra-cell AN leakage Qf NS , where 
Qsns = Qo NS + Qf NS - F° r the SNS AN precoder, Q® NS is obtained as 


qT = Y *»e 

rriy^n 


A A H (h 


\H 


= E 


M M 

tr{ A .„ A "} J 2 P™ = ("t ~ K ) E *», 

m^n rriy^n 

(23) 


where we exploited lf24l Lemma 11] and the independence of A m and li('; m . In contrast, the 
intra-cell AN leakage power is given by CP NS = 


B k E 


h fc A A H (h 


\H 
nnJ 


= B k E 


1 /C A 


H/uk \H 


1 + VrT EA 


= (N T - K)ft 


3 k 

nm 


1 + PtT Em=l A 




as the SNS AN precoder matrix lies in the null space of the estimated channels of all K MTs 
in the n th cell. Similarly, the AN leakage relevant for computation of the POLY data precoder 
is obtained as 

i K 

jfc — • *— *iii-f—iL' (24) 


. K i v^M n 

= (1 - A)Pr Jim 1 £ f3 k PrT Wn : 


- JZ / _ J i nn 

X—>oo A z — J 

k =1 


1 + v T B h 

1 ' Pr 1 Z-jm=l Fn 

2) CNS AN Precoder: For the CNS AN precoder at the n th BS, the AN is designed to lie in 
the null space of the estimated channels between all MK MTs and the BS, i.e., 


k 

nm 


A„ = I 


Nj 1 


H" (H„H"rV, 


(25) 


which has rank L = Nt — MK and exists only if B < 3/M. The corresponding AN leakage to 
the k th MT in the n th cell is given by 


M 


QcNS - ^ Amn E 


m= 1 


A A H (h k ) 

ll mn \ ll mn J 


\H 


M 


1 +PrrEt^n^nl 


(N T -MK)J2ft 


k 

mn 


i+PttE^A 


m= 1 


(26) 


ml 


Furthermore, the CNS AN precoder results in the same Pan as the SNS AN precoder, cf. 
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3) Random AN Precoder: For the random precoder, all elements of A.„ are i.i.d. random 
variables independent of the channel llT3l . i.e., A„ has rank L = N r . Hence, h^ n and A m , Vm, 
are mutually independent, and we obtain 

M r M 


Q random E P T) 


E 


h fc A \ H (h k ) 

ll mn\ ll mn ) 


H 


= Ni 


y B k 

/ j rn 


(27) 


m= 1 


771—1 


Furthermore, we obtain P^ dom = (1 - (j>)P T lim#^ A Ylk=i Pnn- 

Remark 2: If the power and time allocated to channel estimation are very small, i.e., rp T — > 0, 
the S/CNS AN precoders yield the same qQ and Pan as the random AN precoder. This suggests 
that in this regime all considered AN precoders achieve a similar SINR performance for a 
given MT. However, for rp T > 0, the S/CNS AN precoders cause less AN leakage resulting in 
an improved SINR performance compared to the random precoder at the expense of a higher 
complexity. 


B. POLY AN Precoder 


To mitigate the high computational complexity imposed by the matrix inversion required for 
the S/CNS AN precoders, while achieving an improved performance compared to the random AN 
precoder, we propose a POLY AN precoder. Similar to the POLY data precoder, we concentrate 
on the selfish design because of the desired low complexity, and hence, set L = N T — K. The 
proposed POLY AN precoder is given by 

A-n I Nt ^ ^ ^ Vj ^H nn H nn ^ ^ H nre , (28) 

where v — [z/ 0 ,..., uj] T contains the real-valued coefficients of the AN precoder polynomial, 
which have to be optimized. In particular, v> is optimized for minimization of the asymptotic 
average AN leakage caused to all MTs in the n th cell Fan- The corresponding optimization 
problem is formulated as 


min„ Pan = <7 ® 


Tr{G 


a A H a H1 

nn-n-n-^-n ' J nn J 


s.t. :Tr{AfA n } = 1/0-1. 


(29) 


The solution of (l29l ) is provided in the following theorem. 

Theorem 2: For K , N T —> oo, the optimal coefficient vector minimizing the asymptotic average 
AN leakage caused to the users in the n th cell for the AN precoder structure in (1281) is given by 

i/ opt = £ - 1 u;, (30) 
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where [£]jj = ( l+ i +1 + e( l+j and cj = [£ 2 + e(,..., ( J+2 + eC, J+l }. Here, ( 1 denotes again the 

_s_H 

I th order moment of the sum of the eigenvalues of matrix H mt H nn , cf. Theorem 1. e is chosen 
such that Tr{A^A n } = 1//3 — 1. 

Proof: Please refer to Appendix C. ■ 

C. Computational Complexity of AN Precoding 

Similarly to the data precoders, the complexity of the AN precoders is evaluated in terms of the 
number of flops required per coherence interval T. For the SNS AN precoder, the computation 
of A n in (l22l) requires the computation and inversion of a K x K positive definite matrix, 
which entails 0.5(/i 2 + K)(2N T — 1) + K 3 + K 2 + K FLOPs [1281 . and the multiplication 
of an N t x K, an K x K, and an K x N T matrix, which entails N T (N T + K)(2K — 1) 
FLOPs ll28l . Furthermore, the T — r vector-matrix multiplications required for AN precoding 
entail a complexity of (2 N T — 1 )N T FLOPs If28l . respectively. Hence, the overall complexity is 
0.5(A 2 + K)(2N T - 1 ) + K 3 + K 2 + K + N t (N t + K)(2K - 1) + (2 N T - 1 )N t (T - r) FLOPs. 
Similarly, for the CNS AN precoder, we obtain a complexity of 0.5((MA") 2 + MK) (2Nr — 1) + 
(ML) 3 + (MK) 2 + MK + N t (N t + MK) (2MK - 1) + (2 N T - 1 )N t (T - r) FLOPs, whereas 
the random AN precoder entails a complexity of (2 N T — 1 )N t (T — r) FLOPs as only the AN 
vector-matrix multiplications are required. 

Similar to the precoded data vector lt23l Section IV], the POLY precoded AN vector can be 
generated using Horner’s rule. Hence, based on (l28l) . the transmitted AN vector in the n th cell 
can be obtained as 

A n z n z n ^t / oH nn H nri ^z n T —— H nn H. nn (z n T .. .)^ ^ . (31) 

Hence, A n z n can be computed efficiently by first multiplying H nn with z n , which requires 

_£_H 

(2 Nr — 1 )K FLOPs, then multiplying H nn with the resulting vector, which requires ( 2 K — 
1)N t FLOPs, adding z n to the newly resulting vector, and repeating similar operations (J + 1) 
times, see m, mi for details of Homer’s rule. Overall, this leads to a complexity of (J + 
1) ((2 K - l)N T + (2 N t - 1 )K) (T - r) FLOPs. 

V. Comparison of Linear Data and AN Precoders 

In this subsection, we compare the secrecy performances of the considered data and AN 
precoders. Thereby, in order to get tractable results, we focus on the relative performances 
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TABLE I 


SINR OF THE k th MT IN THE n th CELL FOR LINEAR DATA PRECODING AND THE SIMPLIFIED PATH-LOSS MODEL IN (l32l ). 


For this model, Tsrci and Tcrci simplify to Tsrci = F SR ^T, and Tcrci = 




, r crci = 


(1-4>)(3Q+-. 


r SRCI l9 + 1 
0 = 


l+ap T T 


r ORCI 1 ^+ 1 
AND $ = 1+(M-1)PP^ ~ 

l+ap T r 


WHERE 


Data Precoder 

T nk 

SZF 

(l-cfr)0Q+l3cl>(a-6) + (M-l) P 2 8<l>a-P)+P/PT 

SRCI 

1 

l/S(0,P/tfinc,i) + (M-l)p 

CZF 

0<£(l-M/3) 

(1-<A)£Q+/3 < M1-0) + (M-1)p 2 0<M1-M/3)+ j S/P t 

CRCI 

1 

l/G(aP,aP/i'c,nr,l)+(M-l)p 

MF 

9(f) 

(l-d>’)BQ+B<ba+(M-l)p 2 ed>+0/P'r 


of SZF, CZF, and MF fl3l data precoders and SNS, CNS, and random AN precoders. The 
performances of SRCI, CRCI, and POLY data precoders and the POLY AN precoder will be 
investigated via numerical and simulation results in Section EO 

In order to gain some insight for system design and analysis, we adopt a simplified path-loss 
model. In particular, we assume the path losses are given by 


A 


k 

mn 


{ 1, m = n 
p, otherwise 


(32) 


where p e [0,1] denotes the inter-cell interference factor. For this simplified model, a and c in 
([7]) simplify to a = 1 + (M — 1 )p and c = 1 + (M — 1) p 2 . Furthermore, the SINR expressions 
of the linear data precoders considered in Section IIII-AI and the MF precoder considered in lfl3l 
can be simplified considerably and are provided in Table HI where we use the normalized AN 
leakage Q = Q/L. The expressions for the normalized AN leakage 0, the asymptotic average 
AN leakage -Pan* an d the dimensionality L of the considered linear AN precoders are given in 
Table HD 


A. Comparison of SZF, CZF, and MF Data Precoders 

In this subsection, we compare the performances achieved with SZF, CZF, and MF data 
precoders for a given AN precoder, i.e., L and 0 are fixed. Since the upper bound on the 
capacity of the eavesdropper channel is independent of the adopted data precoder, cf. Section 
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TABLE II 

AN LEAKAGE FOR SIMPLIFIED PATH-LOSS MODEL IN lf32l) . 0 AND $ ARE DEFINED IN THE CAPTION OF TABLe[I] 


AN Precoder 

Q 

Pan 

L 

SNS 

{a-6) 

(1 - 4)Fr0 

N t -K 

CNS 

a( 1 — 6) 

(1 — $ 

N t - MK 

Random 

a 

(i - 4>)Pt 

Nt 


III-C1 we compare the considered data precoders based on their SINRs. Exploiting the results in 
Table HI we obtain the following relations between 7 nfc F ’ TnT, and 7 ^ F : 


7nF 

7^| F 


= 1 + P(<rfIT - 1 ) 


and 


7 nfc F 1 a(a l)/3 CZF 

^szf = —-— + —i- —Ink • (33) 


1-/3 1-/3 

Hence, for > yTjf, we require 7 ®f F > 1/c = 1/(1 + p 2 (M - 1)), and for 7 ^ fc ZF > 7 ®f F , we 
need 7 / :zf > l/(pa) = l/[p(l + p(M — 1))]. As expected, (l33l) suggests that for a lightly loaded 
system, i.e., /3 —> 0, all three precoders have a similar performance, i.e., 7 ^ ZF ~ 7 // F ~ 7 /J F . 
In the following, we investigate the impact of the number of MTs and the pilot power on the 
relative performances of the considered data precoders. 

1) Number of MTs: From (l33l) . we find that for 7 ®f F > 7 J)| F and 7 ^ ZF > 7 // F to hold, the 
number of MTs has to meet K < K S z f>mf and K < K c zf>szf, where A" S zf>mf = 

6<j)N T T ^ p(j)6N T 


and K, 


CZF>SZF 


(34) 


(1 - <j>)Q + af + 1/P T ~~ (l~<f>)Q + [a(l-6)+p6M\<i>+l/P T ' 

respectively. Interestingly, both the maximum numbers of MTs for which the SZF data precoder 
is advantageous compared to the MF data precoder, Aszf>mf> and the maximum number of MTs 
for which the CZF data precoder is advantageous compared to the SZF data precoder, Aczf>szf, 
decrease with increasing AN leakage, Q, and increasing number of cells, M, but increase with 
the amount of resources dedicated to channel estimation, p T r (via 9), and consequently with 
the channel estimation quality. However, while Aszf>mf decreases with increasing inter-cell 
interference factor, p (via a), Aczf>szf increases. 

2) Pilot Energy : From (1331) . we find that for 7 // F > 7 /j. F and 7 ^ ZF > 7 ®f F to hold, pilot 
energy p T r has to fulfill 

1 , , 1 


p T r > (Pt't)szf>mf — 


0(l-/3)//3+j 

cl+1/Pt 


— a 


and p T r > (p T r) czf>szf = 


0,+1/Pt 


(35) 


— a 
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where we have assumed that SNS AN precoding is adopted, i.e., Q = a — 9, to arrive at insightful 
expressions. Similar results can be obtained for other AN precoders. From (1351) . we observe that 
MF, SZF, and CZF data precoding are preferable if 0 < p T r < (p r r) szf>mf, (Pt , t)szf>mf < 
p T r < (p t t) czf>szf> and p T r > (jv^)czf>szf, respectively. In general, the more MTs are in 
the system (i.e., the larger /3), the larger the pilot energy has to be to make SZF and CZF data 
precoding beneficial. In fact, from (1151) we observe that if /? exceeds /3 mf = (j)/\a 2 +a/P t+4>— 1], 
MF data precoding is always preferable regardless of the value of p T r. Similarly, if 3 exceeds 
Aszf = Op/ [a 2 + a/P T + 4>p — 1], SZF data precoding is always preferable compared to CZF 
data precoding regardless of the value of p T r. 

B. Comparison of SNS, CNS, and MF AN Precoding 

In this subsection, we analyze the impact of the AN precoders on the secrecy rate. AN 
precoders affect the ergodic capacity of the eavesdropper via L and the achievable rate of the 
MT via the leakage, Q. Since the upper bound on the ergodic secrecy rate of the eavesdropper 
in ([7]) is a decreasing function in L, we have 

zieve | ^ /-yeve | ^ /"yevei ( r ic\ 

^nk I random S O nk |SNS S O n k I CNS- PO) 

On the other hand, from Table UH we observe (Random > Qsns > Qcns- Since according to Table 
H] the SINRs for all data precoders are decreasing functions of Q, for a given data precoder, we 
obtain for the lower bound on the ergodic rate of the k th MT in the n lh cell 

Rnk | random A Rnk | SNS < Rnk I cns- (3V) 

Considering (l36l) . (l37l) . and the expression for the ergodic secrecy rate, Iff; = [R n k — it 

is not a priori clear which AN precoder has the best performance. In fact, our numerical results 

in Section lYll confirm that it depends on the system parameters (e.g. a, (3, M, p T r, and p ) which 

AN precoder is preferable. 

C. Ergodic Secrecy Rate Analysis 

In this subsection, we provide closed-form results for the ergodic secrecy rate for SZF, CZF, 
and MF data precoding for the simplified path-loss model in (l32l) . Thereby, the simplified path- 
loss model is extended also to the eavesdropper, i.e., (3 n E = 1 and /3 m E = />, rn f n, is assumed. 
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Combining ©, ©, and the results in Table HI we obtain the following lower bounds for the 
ergodic secrecy rate of the k th MT in the n th cell: 


Kk > 


log 


log. 


( (Q+l/PT)P+{a-—Q)P<t>+cO<t> ~X<f>+X A 

2 \ (Q+l/Pr)l3+(a-Q)P<l>+(c-l)e<t> (1 ~\)<t>+X ) 

f _(Q+1/ Pr)P+(a-6-Q)p<l>+c0(l-P)<l> ~x<t>+X \ 

2 \ (Q+l/Pr)l3+(a-e-Q)p<l>+(c-l)e(l-P)<t> ' 0-~x)4>+X J 

( (Q+l/PT)P+{a-(i9-Q)P<l>+c6(l-Mp)(t> _ ~X<t>+X 


for MF, 

for SZF, (38) 
for CZF, 


l0g 2 V ( Q+ 1/ /¥)■>'+( u-uH-Q)lH >+0 ~ I- m l-Mt)., ' (1 -x)<t>+x) 

where X = ~ ~ ^,a' r > and Q and L are given in Table [TT] for the considered AN precoders. 
Eq. (l38l) is easy to evaluate and reveals how the ergodic secrecy rate of the three considered data 
precoders depends on the various system parameters. To gain more insight, we determine the 
maximum value of a which admits a non-zero secrecy rate. This value is denoted by a s in the 
following, and can be shown to be a decreasing function of 0 for all conidered data precoders. 
Hence, we find a s by setting = 0 in (l38l) and letting 0 —y 0. This leads to 


a z Q 


a, = 


Qa+c6 Nt / L-\-a/P t 
_(1 -(3)cl 2 0 _ 

Qa-\-cQ(l— (3)Nt / L~\~a/P t 
(1 -M/3)a 2 e 


for MF, 
for SZF, 
for CZF. 


(39) 


w Qa-\-c6(l—M(3 )Nt/L+ a/P t 

Eq. (l39l) reveals that for a given AN precoder, independent of the system parameters, the MF 
data precoder can always tolerate a larger number of eavesdropper antennas than the SZF data 
precoder, which in turn can always tolerate a larger number of eavesdropper antennas than the 
CZF data precoder. This can be explained by the fact that the high AN transmit power required 
to combat a large number of eavesdropper antennas drives the receiver of the desired MT into 
the noise-limited regime, where the MF data precoder has a superior performance compared to 
the S/CZF data precoders. On the other hand, since a s depends on both 0 and L, it is not a 
priori clear which AN precoder can tolerate the largest number of eavesdropper antennas. For 
a lightly loaded network with small 0 and small M, according to Table HU we have L ?« N T 
for all three AN precoders. Hence, in this case, we expect the CNS AN precoder to outperform 
the SNS and random AN precoders as it achieves a smaller Q. On the other hand, for a heavily 
loaded network with large (3 and M, the value of a s of the CNS AN precoder is compromised 
by its small value of L and SNS and even random AN precoders are expected to achieve a larger 


a,. 
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VI. Performance Evaluation 

In this section, we evaluate the performance of the considered secure multi-cell massive MIMO 
system. We consider cellular systems with M = 2 and M = 7 hexagonal cells, respectively, 
and to gain insight for system design, we adopt the simplified path-loss model introduced in 
Section [V] i.e., the severeness of the inter-cell interference is only characterized by the parameter 
p £ (0,1]. The pilot sequence length is r = K. The simulation results for the ergodic secrecy 
rate of the k th MT in the n th cell are based on ©, ©, and the expression for the ergodic rate 
of the MT llT3l Eq. (8)] and are averaged over 5, 000 random channel realizations. Note that, 
in this paper, we consider the ergodic secrecy rate of a certain MT, i.e., the k th MT in the n th 
cell. The cell sum secrecy rate can be obtained by multiplying the secrecy rate of the k th MT 
by the number of MTs, K, as for the considered channel model, all MTs in the n ,h cell achieve 
the same secrecy rate. The values of all relevant system parameters are provided in the captions 
of the figures. To enable a fair comparison, throughout this section, we adopted the selfish SNS 
AN precoder when we compare different data precoders and the selfish ZF data precoder when 
we compare different AN precoders. 

A. Ergodic Capacity of the Eavesdropper for Conventional AN Precoders 

In Fig. U] we show the ergodic capacity of the eavesdropper for the considered conventional 
AN precoders. First, we note that the upper bound in ([7]) is very tight since the number of BS 
antennas is large (Nr = 200) and a < a 2 L / (cNr) holds for all considered AN precoders and 
all consider values of a and (3. Furthermore, as 3 increases, the ergodic capacity of all AN 
precoders decreases since the power allocated to the information-carrying signal of the user that 
the eavesdropper tries to intercept decreases with increasing /3 as the total power allocated to 
the information-carrying signals of all users is fixed. As expected, the eavesdropper’s capacity 
benefits from larger values of a. Furthermore, as predicted in (l36l) . because of their different 
values of L, the CNS AN precoder yields the largest eavesdropper capacity, while the random 
AN precoder yields the lowest. The performance differences between the different AN precoders 
diminish for small values of a and /3 as the dependence of the eavesdropper capacity on L 
becomes negligible for small a, cf. ©, and L « Nr holds for all precoders for small 3, 
cf. Table El 
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3.5 


■ 2.5 


- CNS Upper bound 

O CNS Simulation 
- - SNS Upper bound 
^ SNS Simulation 
- Random Upper bound 
Random Simulation 


Fig. 1. Ergodic capacity of the eavesdropper vs. the normalized number of MTs in the cell, /?, for a system with Nt = 200, 
(j) = 0.75, P T = 10 dB, p = 0.3, and M = 2. 


B. Ergodic Secrecy Rate for Conventional Linear Data Precoders 

In Figs. [2] and [3l we show the ergodic secrecy rates of the k th MT in the n th cell vs. the 
number of BS antennas for the MF, SZF, CZF, SRCI, and CRCI data precoders for a lightly 
loaded and a dense network, respectively, and a fixed power allocation factor of f = 0.75. In both 
figures, the analytical results were obtained from Q, ©, and (fl2l) for the SRCI data precoder, 
(fl6l) for the CRCI data precoder, and (l38l) for the MF, SZF, and CZF data precoders. For all 
considered precoders, the analytical results provide a tight lower bound for the ergodic secrecy 
rates obtained by simulations. Furthermore, as expected, the RCI data precoders outperform the 
ZF data precoders for both the selfish and the collaborative strategies, but the performance gap 
diminishes with increasing number of BS antennas. 

For the lightly loaded network in Fig. [2] we assume M = 2 cells, K = 10 users, and a small 
inter-cell interference factor of p = 0.1. For this scenario, the collaborative designs outperform the 
selfish designs and C/SZF precoding yields a large performance gain compared to MF precoding. 
This is expected from our analysis in Section IV-AI as for the parameters valid for Fig. [2l we 
obtain from (l34l) . Ag ZF>MF ~ 250 and A" CZF>SZF ~ 60 for N T = 400. Intuitively, as the network 
is only lightly loaded, the collaborative data precoder can efficiently reduce interference to the 
other cell despite the pilot contamination. 
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Fig. 2. Analytical and simulation results for the ergodic 
secrecy rate vs. the number of BS antennas, Nt, for a 
lightly loaded network with <j> = 0.75, Pt = 10 dB, 
p T = Pt/K, a = 0.1, K = 10, p = 0.1, and M = 2. 



Fig. 3. Analytical and simulation results for the ergodic 
secrecy rate vs. the number of BS antennas, Nt, for a 
dense network with <j> = 0.75, Pt = 10 dB, p T = Pt/K , 
a = 0.1, K = 20, p — 0.3, and M = 7. 


For the dense network in Fig. [3] we assume M = 7 cells, K = 20 users, and a larger inter-cell 
interference factor of p = 0.3. In this case, for the considered range of N T , the collaborative 
precoder designs are not able to suppress inter-cell interference and AN leakage to other cells 
sufficiently well to outperform the selfish precoder designs. In fact, for Nj = 400, we obtain 
from (l34l) K c zf>szf ~ 16, i.e., our analytical results suggest that the SZF precoder outperforms 
the CZF precoder for K = 20 which is confirmed by Fig. [3} Nevertheless, for Np > 400, the 
ergodic secrecy rate for the CZF data precoder will eventually surpass that for the SZF data 
precoder. 

C. Optimal Power Allocation 

In this subsection, we investigate the dependence of the ergodic secrecy rate on the power 
allocation factor 0 and study the impact of system parameters such as f3, M, and p on the 
optimal 0 that maximizes the ergodic secrecy rate. The results in this subsection were generated 
based on the analytical expressions in ©, ©, and ( ITU ) for the SRCI data precoder, (fl6l) for the 
CRCI data precoder, and (1381) for the MF, SZF, and CZF data precoders. 

Fig. |4] depicts the ergodic secrecy rate of the k th MT in the n th cell for the selfish data 
precoders SRCI, SZF, and MF as a function of the power allocation factor 0. All curves are 
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Fig. 4. Ergodic secrecy rate vs. (f> for different selfish data Fig. 5. Ergodic secrecy rate vs. (j> for different data 

precoders for a network with Pt = 10 dB, Nt = 100, precoders for a network with Pt = 10 dB, Nt = 100, 

p T = Pt/K, a = 0.1, p = 0.1, and M = 7. p T = Pt/ K, a = 0.1, /3 = 0.1, and p = 0.1. 


concave and have a single maximum. For 0 = 0 only AN is transmitted, hence = 0 results 
since no data can be transmitted. For 0 = 1, no AN is transmitted, hence = 0 results 
since the capacity of the eavesdropper becomes unbounded (recall that we make the worst-case 
assumption that the eavesdropper can receive noise-free). For 0 < 0 < 1, a positive secrecy 
rate may result depending on the system parameters and the precoding schemes. Since we keep 
the total transmit power fixed, the transmit power per MT decreases with increasing f3. To 
compensate for this effect, the portion of the total transmit power allocated to data transmission 
should increase. This is confirmed by Fig. [4] where the optimal value of 0 for /3 = 0.5 is larger 
than that for /3 = 0.1. Furthermore, for a given d, the optimal 0 is the larger, the better the 
performance of the adopted data precoder is, i.e., for a more effective data precoder, transmitting 
the data signal with higher power is more beneficial, whereas for a less effective data precoder 
impairing the eavesdropper with a higher AN power is more beneficial. 

In Fig. [5] we show the ergodic secrecy rate vs. 0 for the CRCI, CZF, and SZF precoders. 
Similar to our observations in Fig. |4] for given system parameters, the optimal 0 tends to be 
larger for more effective precoders that achieve a better performance. For the system with M = 7, 
this can be observed by comparing the optimal 0 for the SZF and CZF precoders. Furthermore, 
while for the smaller system with M = 2 cells collaborative precoding is always preferable, for 
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Fig. 6. Ergodic secrecy rate vs. tj> for different AN Fig. 7. a a vs. p for different data and AN precoders for 

precoders for a network with Pt = 10 dB, Nt = 100, a network with Pt = 10 dB, Nt = 100, p T = Pt/K , 

p T = Pt/K, M = 2, p = 0.1, and a = 0.1. p = 0.3, and M = 2. 


M = 7, SZF precoding outperforms CZF and CRCI precoding for all considered values of (f>, as 
the collaborative designs are not able to effectively suppress the interference and AN leakage to 
the (M — 1 )K = 60 users in the other cells with the available Np = 100 antennas. In particular, 
from (l34l) . we obtain A" czf >szf < 18 for M = 2 and Aczf>szf < 5 for M = 7, which confirms 
the results shown in Fig. [5] Fig. |6] depicts the ergodic secrecy rate vs. o for the considered 
conventional AN precoder structures. We consider a lightly loaded network with /3 — 0.2 and a 
moderately loaded network with (3 = 0.4. For 3 = 0.2, the CNS AN precoder outperforms the 
SNS AN precoder since, in this case, for the CNS AN precoder, the negative impact of having 
(slightly) fewer dimensions available for degrading the eavesdropper’s channel (smaller value of 
L) is outweighed by the positive impact of causing less AN leakage (smaller value of Q). On 
the other hand, for /3 = 0.4, the CNS AN precoder has a substantially smaller L than the SNS 
precoder which cannot be compensated by its larger Q. Despite having the largest value of L, 
the random AN precoder has the worst performance for both considered cases because of its 
large AN leakage. 
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D. Conditions for Non-Zero Secrecy Rate 

In Section IV-Cl we showed that a positive ergodic secrecy rate is possible only if o < a s . In 
Fig. [71 using (l39l) . we plot a s as a function of /3. In the left hand side subfigure, we compare 
MF, SZF, and CZF data precoding for SNS AN precoding, and in the right hand side subfigure, 
we compare random, SNS, and CNS AN precoding for SZF data precoding. The comparison 
of the data precoders reveals that although SZF and CZF entail a much higher complexity, MF 
precoding achieves a larger a s . Therefore, if the eavesdropper has a large number of antennas 
and small ergodic secrecy rates are targeted, simple MF precoding is always preferable. On the 
other hand, whether SNS or CNS AN precoder is preferable depends on the system load. For 
small values of /3, CNS AN precoding can tolerate more eavesdropper antennas, whereas for 
large values of /3, SNS AN precoding is preferable. Random AN precoding is outperformed 
by SNS AN precoding for any value of j3. A closer examination of (l39l) reveals that this is 
always true if S/CZF data precoders are employed. However, for the MF data precoder, there 
are parameter combination for which random AN precoding outperforms SNS and CNS AN 
precoding. 

E. Low-Complexity POLY Data and AN Precoders 

In this subsection, we evaluate the ergodic secrecy rates of the proposed low-complexity 
POLY data and AN precoders. To this end, we consider again a lightly loaded network with 
little inter-cell interference (M = 2, 13 = 0.1, p = 0.1) and a dense network with more inter-cell 
interference (M = 7, f3 = 0.15, p = 0.3). All results shown in this section were obtained by 
simulation. For each simulation point, the optimal value of <fi was found numerically and applied. 
In Figs. [8] and [9[ we show the ergodic secrecy rate of the k th MT in the n lh cell as a function 
of the pilot energy, rp T . As expected, for all considered schemes, the ergodic secrecy rate is 
monotonically increasing in the pilot energy since more accurate channel estimates improve the 
performance. 

In Fig. [8} we depict the ergodic secrecy rates for the proposed POLY data precoder for 
different values of X and compare them to those of conventional selfish data precoders. For 
the sake of comparison, all data precoders are combined with the SNS AN precoder. As the 
number of terms of the polynomial X increase, the performance of the POLY data precoder 
quickly improves and approaches that of the SRCI data precoder. The convergence is faster for 
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M = 2,p = 0.1,/3 = 0.1 M = T,p = 0.3,/5 = 0.15 



Fig. 8. Ergodic secrecy rate for POLY and conventional 
selfish data precoders for a network employing the optimal 
<f>, Pt = 10 dB, Nt = 200, and a = 0.1. 


M = 2,p = 0.1, (3 = 0.1 M = 7,p = 0.3,/3 = 0.15 



Fig. 9. Ergodic secrecy rate for POLY, SNS, and random 
AN precoders for a network employing the optimal <p , 
P T = 10 dB, Nt = 200, and a = 0.1. 


the dense network considered in the right hand side subfigure, where the performance difference 
between all precoders is smaller in general since interference cannot be as efficiently avoided 
as for the lightly loaded network. 

In Fig. |9l we show the ergodic secrecy rates for the proposed POLY AN precoder for 
different values of J and compare them to those of the random and SNS AN precoders. For 
the sake of comparison, all AN precoders are combined with SZF data precoding. The POLY 
AN precoder quickly approaches the performance of the SNS AN precoder as the polynomial 
order J increases. Similar to the POLY data precoders, the convergence is faster for the dense 
network where the performance differences between different AN precoders are also smaller. 
For the denser network, even the random AN precoder is a viable option and suffers only from 
a small loss in performance compared to the SNS AN precoder. 

F. Complexity-Performance Tradeoff 

In this subsection, we investigate the tradeoff between the ergodic secrecy rate performance 
and the computational complexity of the proposed data and AN precoders in Figs. OH and ITT! 
respectively. In particular, Figs. [TO] and QT| depict the ergodic secrecy rate on the left hand side 
and the computational complexity (in Giga FLOP) on the right hand side, both as a function of 
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Ergodic secrecy rate Computational complexity 



Fig. 10. Ergodic secrecy rate (left hand side) and 
computational complexity (right hand side) of various 
linear data precoders for a network employing Pt = 10 
dB, N t - 1000, p T = Pt/K, M = 2, p = 0.1, 
T — t = 100, and an SNS AN precoder. 


Ergodic secrecy rate Computational complexity 



Fig. 11. Ergodic secrecy rate (left hand side) and 
computational complexity (right hand side) of various 
linear AN precoders for a network employing Pt = 10 
dB, N t = 1000, p T = P T /K, M = 2, p = 0.1, 
T — t = 100, and an SZF data precoder. 


the numbers of users in a cell. For the considered setting, the performance gains of collaborative 
data and AN precoding compared to selfish strategies are moderate, but the associated increase 
in complexity is substantial, especially for large K. 

Fig. ITOl illustrates that for the considered setting a POLY data precoder with X = 1 achieves a 
better performance than the MF precoder but has substantially lower complexity than the SRCI 
precoder. For large X, the POLY data precoder has a lower complexity than the SRCI precoder 
for large K. However, even for small K , the POLY precoder may be preferable as it does not 
incur the stability issues that may arise in the implementation of the large-scale matrix inversions 
required for the SRCI precoder. 

Fig. QU shows that for the considered setting the proposed POLY AN precoder with J = 1 
outperforms the random AN precoder. The POLY AN precoder with J — 5 achieves almost 
the same performance as the SNS AN precoder but with a substantially lower complexity. We 
further observe that for small K , because of its efficient implementation via Horner’s scheme, 
cf. dm the proposed POLY AN precoder requires an even lower complexity than the random 
AN precoder. 
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VII. Conclusion 

In this paper, we considered downlink multi-cell massive MIMO systems employing linear 
data and AN precoding for physical layer security provisioning. We analyzed and compared the 
achievable ergodic secrecy rate of various conventional data and AN precoders in the presence of 
pilot contamination. To this end, we also optimized the regularization constants of the selfish and 
collaborative RCI precoders in the presence of AN and multi-cell interference. In addition, we 
derived linear POLY data and AN precoders which offer a good compromise between complexity 
and performance in massive MIMO systems. Interesting findings of this paper include: 1) 
Collaborative data precoders outperform selfish designs only in lightly loaded systems where 
a sufficient number of degrees of freedom for suppressing inter-cell interference and sufficient 
resources for training are available. 2) Similarly, CNS AN precoding is preferable over SNS AN 
precoding in lightly loaded systems as it causes less AN leakage to the information-carrying 
signal, whereas in more heavily loaded systems, CNS AN precoding does not have sufficient 
degrees of freedom for effectively degrading the eavesdropper channel and SNS AN precoding 
is preferable. 3) For a large number of eavesdropper antennas, where only small positive secrecy 
rates are achievable, MF data precoding is always preferable compared to SZF and CZF data 
precoding. 4) The proposed POLY data and AN precoders approach the performances of the 
SRCI data and SNS AN precoders with only a few terms in the respective matrix polynomials 
and are attractive options for practical implementation. 

Appendix 

A. Proof of Proposition 1 

Considering © and ©. the effective signal power, i.e., the numerator in ©, can be expressed 
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where L„ >fc = (H n „H" - (h£ n ) H h£ n + ^1^) \ X nk = E[h k nn L n , k (h^ n ) H ], and A nk = 

E[h^ n L n fc (h^ n )^]. On the other hand, the intra-cell interference term in the denominator of 
© can be expressed as 
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where H n>fc is equal to U nn with the k th row removed, and Y nk = E[h^ n L rufc H r ^H nJ X n:fc (h^ r J 
and B nk = E[h^ m L n>fc Hj fc H n)fc L n)fc (h^ n ) ff ]. 

Due to pilot contamination, the data precoding matrix of the m th BS is a function of the channel 
vectors between the m th BS and the k th MTs in all cells. Hence, the inter-cell interference from 
the BSs in adjacent cells is obtained as 


E[|hLU 2 ] = 


li(X nk + A nt ) 2 1 +PrTY. 


M 


ok 
l^m oml 


(42) 


(1 + W*) 2 1 + r ,rE,'=,L 

Meanwhile, by exploiting (l40l) . (l42l) . and the definition of the variance, i.e., var [x = E[x 2 ] — 
E 2 [x], we obtain for the first term of the denominator of ©, var[hjj n f nfc ] = 

1 +PtT z^m=l Pnm 


According to |20l Eq. (16)] and [27] Theorem 7], for N r —y oo and constant /3, X nk converges 
to G((3, h i) defined in (fill) and A nk —>■ 0. Similarly, Y nk and B nk approach 
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respectively, where i) = 

Moreover, the inter-cell interference from other MTs (i.e., not the k th MTs) is calculated as 


E 


h m nFm,fcF mjfc (h mrl ) 


H 


= E 


t r 


= K - 1, 


(45) 


where F m is equal to F m with the k th column removed. The first equality in (|45T) is due to 
the fact that the precoding matrix for the other MTs (i.e., not the k th MTs) in adjacent cells are 
independent of h] m and 11241 Lemma 11], while the second equality holds for N T —y oo. 

On the other hand, the constant scaling factor ^i for SRCI precoding is given by 11201 Eq. 

( 22 )] 

(46) 


7i = 


GWi Hi) + Hi-q^-G{( 3, Hi) 


Hence, employing (I40l)-(l46l) in ©, the received SINR in (fTOl) is obtained, which completes the 
proof of Proposition 1. 
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B. Proof of Theorem 1 

The objective function in (l20l) can be rewritten as 

r x /„ * h \ i+1 x 
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1 + S -Pan + cTr {S n } , (47) 


where we exploited E[s „s?] = I K , the definition of Pan given in Theorem 1, the definition of 
Fn in (1 1 81), the definition ^^ U ^ firm T Hnra^ and H nri Nt * 

In the following, we simplify the right hand side (RHS) of (1471) term by term. To this end, 
we denote the first three terms on the RHS of (l47l) by t\, t 2 , and t 3 , respectively. Using a result 
from free probability theory ll26l . the first term converges to ll23l Theorem 1] 
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as matrix D nn is free from Jf i=0 //, ( H nn H nn ) . Similarly, the third term converges to 
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Furthermore, the second term can be rewritten as 
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where (a) follows again from [[23] Theorem 1] and (b) results from ETr {H nn D n ,,H, m }] = 
Tr {D rm A„}, where A n is defined in Theorem 1, (fl8l) . and the constraint in ([20b . 

_~_H 

Exploiting (l48l)-(l50l) and the eigen-decomposition of matrix H n?l H ren = TAT H , where diag¬ 
onal matrix A = diag (Ai,..., A+ contains all eigenvalues and unitary matrix T contains the 
corresponding eigenvectors, the asymptotic average MSE becomes 
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Next, we introduce the Vandermonde matrix C\ E R Kx( - I+1 \ where |C | \ hJ = 1 , and A = 

[Ai,..., A k] t , which allows us to rewrite (1511) in compact form as 


mse n = lim — E 
K —}oo K 


<? 2 pTr {D nn } n T C? A 2 C lM - 2<,VpTr {D^ 2 } /x T Cf A 
+1 + g 2 P A N + ? 2 Tr {5] n } + (fpNxTi {D nn A„} . 

Similarly, the constraint in (l20l) can be expressed as 
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(53) 


Thus, the Lagrangian function of primal problem (l20l) can be expressed as £, (ji. q) = mse n + 
e^lim^oo ^E[/it T Cf ACi/x] — Nt), where ei is the Lagrangian multiplier. Taking the gradient 
of the Lagrangian function with respect to /r, and setting the result to zero, we obtain for the 
optimal coefficient vector /x opt : 
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Furthermore, taking the derivative of £i(/z, q) with respect to q and equating it to zero, and 
multiplying both sides of (l54l) by /x T and applying (1531) . we obtain 

(55) 
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The expressions involving Ci, A, and A in (1541) can be further simplified. For example, we obtain 
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(l54l) in a similar manner anc 
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. Simplifying the other terms in 


inserting (1551) into (l54l) we obtain the result in Theorem 1. 


C. Proof of Theorem 2 

Exploiting E[z n z^] = I Nt , the constraint in (l29l) . and a similar approach as was used to arrive 
at (l24b . the objective function in (|29|) can be simplified as P A n = 
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Using (1281) and a similar approach as in Appendix B, (1561) can be rewritten as 
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Defining Vandermode matrix C 2 G E^ x (- /+| ), where [C 2 ]ij = A { 1 , we can rewrite (1571) in 
compact form as P AN = 


qN T Tr {D nn } lim —E 

K—too K 


- 2v T &> A A + 1 T A + u T Ci A 3 C 2 j/ 


+ (1 - 4>)P t Tr{D nn A n }, (58) 


where 1 denotes the all-ones column vector. Taking into account the constraint in (l29l) . we can 
formulate the Lagrangian as C, 2 {y) = P AN + 6 2 (lim.A'—Kxj ^E^C^A 2 ^ - 2 I / T C[A] + 1) 
with Lagrangian multiplier e 2 . The optimal coefficient vector i/ opt is then obtained by taking the 
gradient of the Lagrangian function with respect to v and setting it to zero: 
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K—t oo 


C[A 2 (A + elx) C 2 


v = 
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K —>-00 


C£(A + eI*)A , 


(59) 


where we used e = )/; . VtT ^| D — y. Simplifying the terms in (l59l) by exploiting a similar approach 
as in Appendix B, we obtain the result in Theorem 2. 
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